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Abstract 



We present a new approximation scheme that aUows us to increase 
the accuracy of analytical predictions of the power spectra of inflation- 
Q ' ary perturbations for two specific classes of inflationary models. Among 

^ I these models are chaotic inflation with a monomial potential, power-law 

C/3 . inflation and natural inflation (inflation at a maximum). After reviewing 

^ ' the established first order results we calculate the amplitudes and spec- 

tral indices for these classes of models at higher orders in the slow-roll 
. , - parameters for scalar and tensor perturbations. 

>< 

^ • 1 Introduction 

Inflationary cosmology Q is facing exciting times due to a new generation of 
ground and satellite based experiments to be carried out (e.g., the SDSS, MAP 
and Planck experiments (2|). Upcoming observations will allow us to determine 
the values of cosmological parameters with high confidence. To measure the 
values of these parameters it is necessary to make assumptions on the initial 
conditions of the density fluctuations that evolved into the observed large scale 
structure and CMB anisotropics. 
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In the simplest inflationary scenario a single scalar fleld (p drives the accel- 
erated expansion of the Universe. Cosmological perturbations are generated by 
quantum fluctuations of this scalar fleld and of space-time. Perturbations are 
adiabatic and gaussian and are characterized by their power spectra. Usually 
these spectra are described in terms of an amplitude at a pivot scale and the 
spectral index at this scale. For general models, these quantities are difficult 
to compute exactly. The state-of-the-art in such calculations are the approxi- 
mated expressions due to Stewart and Lyth j^, which are obtained up to next- 
to-leading order in terms of an expansion of the so-called slow-roll parameters. 
This expansion allows to approximate the solutions to the equations of motion 
by means of Bessel functions. A recent analysis of Wang et al. [Q shows that 
the Bessel function approximation cannot be further improved. We will show 
below that this result depends critically on the assumptions that are made on 
the relative order of the various slow-roll parameters and that there are two 
regions in the slow-roll parameter space, where higher accuracy results may be 
obtained. 

To reliably compare analytical predictions with measurements, an error in 
the theoretical calculations of some percents below the threshold confldence 
of observations is required. In Ref. |^ it was shown that amplitudes of next- 
to-leading order power spectra can match the current level of observational 
precision. However, the error in the spectral index and the resulting net error in 
the multipole moments of the cosmic microwave background anisotropics might 
be large due to a long lever arm for wave numbers far away from the pivot scale 
A clever choice of the pivot point is essential [|[ |^ for todays and future 
precision measurements. The slow-roll expressions as calculated in ||] are not 
precise enough for the Planck experiment Recently, Stewart and Gong Q 
presented a new method, also based on the slow-roll expansion, that allows to 
obtain analytical expressions at higher orders in the slow-roll parameters. 

In this paper we show that the Bessel function approximation can be im- 
proved to higher orders, without being in conflict with the general argument of 
[H . We calculate the amplitudes and indices of the scalar and tensorial pertur- 
bations up to higher orders for two speciflc classes of inflationary models. The 
first class contains all models that are 'close' to power-law inflation, one example 
is chaotic inflation with y c>c Our second class of models is characterized 
by extremely slow rolling, an example is inflation near a maximum. 

2 The standard formulas 

Let us quickly review the derivation and results of the Bessel function approxi- 
mation. The slow-roll parameters are defined as, [^, 



<0) ^ - 



K H 



2 IH' 




(1) 
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with the equations of motion 



|: = 2e(e-ry), ^=eil-e- (2) 

H is the Hubble rate, dot and prime denote derivatives with respect to cosmic 
time and (j), n = 87r/mpj, and mpi is the Planck mass. By definition e > and 
it has to be less than unity for inflation to proceed. 



2.1 Amplitudes of inflationary perturbations 

We call standard those formulas which are considered to be the state-of-the- 
art in the analytical calculation of perturbations spectra, i.e, those obtained by 
Stewart and Lyth 

The general expression for the spectrum of the curvature perturbations is P| 



Ufc(r) are solutions of the mode equation M, nl 



dr^ \ z dr 

where r is the conformal time and z is defined as z = a(j)/H. The potential of 
the mode equation (^) reads ||, |^ 

Despite its appearance as an expansion in slow-roll parameters, Eq. (|^) is an 
exact expression. 

The crucial point in the Stewart and Lyth calculations is to use the solution 
for power-law inflation (where the slow-roll parameters are constant and equal 
each other) as a pivot expression to look for a general solution in terms of 
a slow-roll expansion. An answer to whether the slow-roll parameters can be 
regarded as constants that differ from each other is found it looking at the exact 
equations of motion (|^) . 

The standard procedure is to consider the slow-roll parameters so small that 
second order terms in any expression can be neglected. Now aH in (|^) can be 
replaced with help of {aH)~^ ~ "'''(1 ~ e)- Then, according with Eq. (Q) the 
slow-roll parameters can be fairly regarded as constants and Eq. (|^) becomes 
a Bessel equation readily solved. From that solution the scalar amplitudes are 
written as 



nir ' m% \H'\ 



(6) 
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where v is given by 



1 - e 



(7) 



and fc is the wavenumber corresponding to the scale matching the Hubble radius. 
Expanding solution with v given by z/ = 3/2 + 2e — 77 and truncating the 
results to first order in e and r\ the standard general expression for the scalar 
spectrum is obtained, 



(8) 



where C = —2 + ln2 + 7 ~ —0.73 is a numerical constant, and 7 is the Euler 
constant that arises when expanding the Gamma function. Eq. (j^) is called 
the next-to-leading order expression for the spectrum amplitudes of scalar per- 
turbations, and from it the leading order is recovered by neglecting first order 
terms for e and r\. 

The corresponding equation of motion for the tensorial modes is 



where, 



10 = 2a^i/^fl-i. 



(9) 



(10) 



Neglecting any order of e higher than the first one, Eq. can be written as 

^ ' " - (11) 



1 d'^a 
a dr 



2 ^2 I 



where. 



u = h - ~ - + e . (12) 

^ 1 - e 2 2 ^ ^ 

This way, Eq. (|) can be also approximated as a Bessel equation with solution. 



12k 



-H 



(13) 



k=aH 



Substituting /i given by Eq. (|12[) in Eq. (|13|), expanding on e and truncating to 
first order, the Stewart-Lyth next-to-leading order result is obtained 



2k 



(14) 



k=aH 



The leading order equation is recovered by neglecting e. 
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2.2 The spectral indices 

To derive the expressions for the spectral indices we introduce here a variation 
of the standard procedure that has the advantage to allow a careful bookkeeping 
of the order of slow-roll expressions. 

First, let us assume the following ansatz for the power spectrum of general 
inflationary models, 

7^^1/2 ^ Me, r;,e)/(e,ry,0, (15) 

where /(e, 77, ^) may be written as a Taylor expansion while h{e, rj, ^) is a general 
function which may be singular or discontinuous at e = 0. Note that any 
function can be decomposed into this form. Wc must also consider the following 
function of the slow-roll parameters: 

K H d(f> 1 _ 

In this expression all of the parameters are to be evaluated at values of cor- 
responding to fc = aH. The functions / and g can be expanded in Taylor 
series, 

/(e,77,^) = floo + flioe + 020?? + 030^ 

+ fliie^ + 012677 + 0136,^ + 022?]^ + 02377.^ + 033^^ ••• , (17) 
g{e,7j,0 - -(l + e + e^-f ...). (18) 

We proceed with the derivation of the equations for the scalar spectral index 



From Eq. ( [l5| ) we obtain 
dlnVR^^'^ _ rfln/i(e,77,C) d 



2 3 
where 

-P(e-'7:C) = aioe+a20'7 + a3oC+aiie^+ai2e??+ai3e^+a22?7^+«23??C+a33C^H , 

and (iij = aij/aoQ. After differentiation, 

^^^^if- = ^ ^ + - ^' ^) + ^' ^) - • • •) ^' ' 

(20) 

where, 

P' = aioe' -1-02077' -Fa3oC'+2aiiee'-|-ai2 (e7;)'-fai3 (eO'+2a227/77' -1-023 ('7C)'+2a33^C' 
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plus higher order derivatives. Using Eqs. ( p^ ) and (|T^), and the definitions of 
the slow-roll parameters (□) we obtain that 



ns -1 



24 + (i-P + P^-^^.)/- 



{l + e + e'^ + + ■■■), (21) 



where P' is written now as, 

V K 

(22) 

Expression ( pi[ ) can be used to any order whenever information on the function 
h and the coefficients of the expansion ( p7| ) is available. The standard result 
of Stewart and Lyth, Eq. (|^), has been tested as a reliable approximation for 
the scalar spectrum of general infiationary models |^ . Then, it is reasonable to 
assume 



1, aio--(2C+l), 



1 U H 

2^V 27e ' 
C. 



(23) 
(24) 



With these assumptions, and ignoring the assumptions (constant slow-roll pa- 
rameters) used to obtain Eq. (ph, it is obtained that 



2 ^h' 

-yfe— = -2e + f]., 

K ll 



(25) 



and to second order (the order for which information on the coefficients for (g2| 
is available from Eq. (H)), expression ( pT| ) reduces to 



71S - 1 = -4e + 2r; - 8(C + l)e^ + 2(5C + 3)e?7 - 2C^^ 



(26) 



which is the standard result for the scalar spectral index. 

The equation of the tensorial index can be derived in a similar way. By 
definition 



dlnVg 



(27) 



The standard equation for tensorial modes is given by Eq. (p^). If we now use 
expression ( pl| ) substituting ng — 1 by nx and 



'2k 



then we obtain 



aio 



2 ^h' 

K h 



-ic + i), 



and the standard expression, 

nr = -2e - 2(2C + 3)6^ + 4(C + l)e?/ . 



(28) 
(29) 

(30) 
(31) 
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3 Generalizing the Bessel approximation 



As can be observed from Eq. (|^) , the standard assumption used to approximate 
the slow-roll parameters as constants can not be used beyond the linear term 
in the slow-roll expansion. Hence, from this point of view, the feasibility of 
using the Bessel equation to calculate the power spectra is limited to this order. 
Nevertheless, what it is actually needed is the right hand sides of Eqs. (|^) being 
negligible. That can also be achieved if |u| = |e — 77I ^ e (we shall call this 
the power-law approximation) or if e ^ |ti| (we call this the extreme slow- 
roll approximation). In both cases we require luff^^l = — -I- 3eu| <JC 
min(e, \v\). No discrepancy arises with the conclusions in Ref. Q], where 

dln(e) 
dN 

is assumed, N being the number of e-foldings. As explained above, this condition 
is not fulfilled in neither of our approximations. 



3.1 Power-law approximation 

Power-law inflation is the model which gives rise to the commonly used power- 
law shape of the primordial spectra, although not always properly implemented 
[0 . The assumption of a power-law shape of the spectrum has been successful in 
describing large scale structure from the scales probed by the cosmic microwave 
background to the scales probed by redshift surveys. It is reasonable to expect 
that the actual model behind the inflationary perturbations has a strong simi- 
larity with power-law inflation. For this class of potentials the precision of the 
power spectra calculation can be increased while still using the Bessel approxi- 
mation. If we have a model with w oc e", then, according with Eq. (0), the first 
slow-roll parameter can be considered as constant if terms like e"^^ and with 
higher orders are neglected. These considerations, plus the condition ^ 
imply the right hand side of the second equation in (||) to be also negligible 
and, this way, 77 can be regarded as a constant too. The higher the order in the 
slow-roll parameters, the smaller should be the difference between them, finally 
leading, for infinite order in the parameters, to the case of power-law inflation. 
That means that this approach to the problem of calculating the spectra for 
more general inflationary models is in fact an expansion around the power-law 
solution. An example is given by chaotic inflation with the potential V oa (f)" 
with a > 2. In this case the slow- roll parameters are given by e ~ a/ (4iV) and 
V ~ 1/(2A^) [0, i.e., e> V. N denotes the number of e-folds of inflation. 

Let us proceed with the calculations. In fact all that has to be done is to 
repeat the calculations of the previous sections keeping the desired order in the 
slow-roll expansions and neglecting all the terms similar to the right hand sides 
of Eqs. (^). For example, the next-to-next-to-leading order expression for the 
scalar power spectrum will be obtained from solution (|6|) but now with v given 
by J/ = 3/2 + e + t; + e^. Expanding and keeping terms up to second order, the 
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final expression is 



(32) 



k=aH 



where B = -2 + + 7r^/4 ~ 0.73. Tliis result is consistent with the result 
obtained in Ref . [|| . Eq. ^ can be readily recovered from Eq. (|3^) by neglecting 
second order terms. For the corresponding expression of the tensorial power 
spectrum, /i must be = 3/2 + e + and, 



=[l-iC + l)e +{B- l)e2] 



'2k 



-H 



(33) 



k=aH 



Equation is obtained from Eq. (|3^) by neglecting second order terms of e. 



3.2 Extreme slow-roll approximation 

There is no reason to reject the possibility of an inflationary model with \v\ 3> e. 
As we shall see later, some important models belong to this class. For these 
cases, the accuracy of the calculations of the amplitudes can also be increased. 
We shall focus on the next-to-next-to-leading order. With regards of conditions 
(||), we neglect terms like e^, ev and — but we keep terms. Repeating 
the calculations under this set of assumptions we obtain for the amplitudes of 
the scalar spectrum: 



V,'/^^^[l~iC+l)e-Cv + Bv ^ 



. (34) 

H 



which again agrees with the result from g in the appropriate limit. For tensorial 
perturbations we note that the Bessel function index ( |l2|) does not depend on 
V hence, no term like v'^ will arise at any moment of the calculations. This way, 
the expression for the tensorial amplitudes is given by Eq. (^. 



3.3 The spectral indices 

In general, if the approximation neglecting the right hand sides of Eqs. 
is consistently taken into account, expressions for the spectral indices at any 
order n are easy to be derived noting that in these cases expression ( p2| ) always 
vanishes and, 

(1 + e + e^ + ... + 6-^), (35) 

where, rii is ns — 1 or tit and hi is correspondingly given by Eqs.(^3|) and (p8|). 
This way it is obtained 

ns{k)-l = (-4e + 277)(l + e + e2 + ... + e") 

= -2v - 2e(l + e + e^ + ... + e") , (36) 
nrik) = -2e(l + e + £2 + ... + £»). (37) 



K hi 
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3.3.1 Relaxing the assumptions 

In the same way as the spectral indices ( p6| ) and (|3|) have been derived up 
to second order in the slow-roll parameters, we can now relax our assumptions 
on the slow-roll parameters and include higher orders beyond those already 
included in (^) and (37). 

Let us first consider the case of the power-law approximation. The ampli- 
tudes are given by Eqs. ( |3^ ) and (p3|). We now drop the assumption that e,f, 
and are constant and keep all the contributions from the derivatives of these 
terms, i.e., ev^v/H^e^v. However, we have to make sure that no derivatives 
of terms that have been neglected in the amplitude show up, thus we have to 
neglect the derivatives of eu,u^, and v/H, which gives rise to the conditions 
et;2 = -evl{2H),vv/H = 0, and [A{v / H) / At] / H = 0. This means that we are 
allowed to keep the terms e"^, e^v, and ev/H beyond the standard expression for 
the spectral indices. 

Consistently with our approach, we must write the expressions for P and P' 
in terms of our basic parameters e and v. These expressions are, 

P = aioe + a20t^ + aiie^ + ai2eu -I- a22U^ + ■ ■ ■ , (38) 
P' = 2aiQev + a2nv / H 

+ Ahiie^v + ai2e(2w2 + v/H) + 2a22Vv/H +■■■ , (39) 



where y^2e/ra/ — v/H — 2e^ — 'ierj + For the scalar contribution in the 
power-law approximation 

5io = -(C+l), a2o = -C, dn^B-1. (40) 

Using Eqs. (|3^) and (|3^), and taking into account the relaxed approximations 
discussed above, expression ( ^T| ) reduces to 

ns-l = ~2e-2v-2€^ -2{2C + 3)€v-2Cv/H 

-2e^ -2{6C +17 -TT^)e^v -2Cev/H . (41) 

Correspondingly, for the tensorial index 

aio^-{C+l), (42) 

and 

riT = -2e - 2e2 _ 4(C + l)eu -26^-2(6(7+ 16 -7r2)e2u. (43) 

In a similar manner we can derive spectral indices for the extreme slow-roll 
approximation. Now we have to make use of the conditions e^v — 0,ev/H = 
—2ev^, and [d{v/ H)/dt]/ H = 0. We have now for the scalars 

5io = -(C + l), = -C, 522 = 5, (44) 

and the spectral index reads, 

ns-l = -2e-2v-2e^ -2{2C + 3)ev-2Cv/H 

-2e' + ACev"^ - (8 - TT^)vv/H , (45) 
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while the tensorial index is identical to the standard second order slow-roll 
result, 

UT = -2e - - 4(C + l)ei; . (46) 

To conclude this section let us not that our expressions for the spectral 
indices are in full agreement with the results of Stewart and Gong |8|, taking 
the corresponding approximations into account. This is a quite nontrivial test 
of our results and of the results of Ref . |^ . 



4 Testing the expressions 

To test all the expressions presented in this paper, we will use the well known 
exact results for amplitudes and indices calculations in the cases of power-law 
and natural inflation. 

Power-law inflation [ pT| is an inflationary scenario where, 



(47) 



a(t) (X , H(<p) cx exp — . / — ( 

' " 2p 



y (0) cx exp — W — ( 
p 



with p being a positive constant. It follows from (|l|) that in this case the slow-roll 
parameters are constant and equal each other. 



v = ^ = ^/p- 



(48) 



Thus, this is the limit case of the power-law approximation. Note that condition 
e < 1 implies p > I. 

For this model, the power spectrum of scalar perturbations is given by p2[ 



1-1 
p 



\m 



k=aH 



with 



1 



1 



Expanding up to second order for p ^ 1 it is obtained that 



K 

47r 



l_(C-f 1)1-^(5-1)1 
p 



(49) 



(50) 



in full correspondence with Eq. ( p^ ) when relations (^8|) are taken into account. 
Testing the tensorial amplitudes can be easily done by checking that Eqs. ( ^2| ) 
and (^3|) indeed satisfy the relation between amplitudes characteristic of power- 
law inflation, i.e.. 



Vp 



(51) 
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For the spectral indices one can see that substitution of relations i 
and @, as weU as in Eqs. Q and (Esl), yield 



riQ — 1 jiT if 1 1 



P 



p 



1 



1-p 



j) in Eqs. (g6|) 
(52) 



the corresponding relation between indices in power-law inflation. 

Another inflationary scenario where precise predictions for the spectra am- 
plitudes and indices can be done is natural inflation p^ . In this case, the 
potential is given by, 

1/ = A4 1±cos(|) , (53) 

where A and / are mass scales. For simplicity we choose the plus sign in the 
remaining calculations. The point here is to analyze inflation near the origin 
so that the small-angle approximation applies, i.e., </) <C /. In this case the 
slow-roll parameters are. 



1 



2 1 
3^ 



1 + ^ 



2 1 
3^ 



(54) 

(55) 

(56) 
(57) 



As it can be observed, in this approximation natural inflation belongs to the 
class of models that fulfill the conditions of the extreme slow-roll approximation. 
With the values of the parameters given above, Eq. (^) can be used with the 
corresponding ly = 3/2 + v. After expanding and truncating at the proper order 



is obtained, consistent with Eq. (|3j 
Eq. (^ and (^) are reduced to 



(58) 



k=aH 

With this same set of assumptions. 



ns - I = 21], 
the well known result for natural inflation. 



(59) 



5 Conclusions 

In this paper we introduced a simple procedure which allows to improve the 
precision of the predictions of inflationary perturbations for two specific classes 
of inflationary models. The method is based on neglecting the evolution in 
time of the slow-roll parameters e and v = e — r/. All monomials of the slow- 
roll parameters e and v that are larger than 'niax{e / H , v / H) can be taken into 
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account, the rest is dropped consistently. Two cases arise. In the first one, e 
is greater than v. We called this approach power-law approximation. In the 
opposite case, the value of e is smaller than v (extreme slow- roll approximation). 
These approximations plus the standard approximation cover a large space of 
inflationary models. 

We briefly discussed the standard approach to the derivation of the spectral 
indices. It was noted that the information contained in the expressions for 
the power spectra amplitudes obtained within the Bessel formalism is used as 
an approximation to derive the spectral indices for more general models. The 
obtained expressions have been tested against the exact results for power-law 
and natural inflation. 
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